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ABSTRACT 

We show that simple kinematic arguments can give limits on the timescale 
and amplitude of variabilities in gamma-ray burst (GRB) afterglows, especially 
when the variability timescale is shorter than the observed time since the burst 
At < t. These limits help us to identify the sources of afterglow variability. The 
afterglows of GRB 011211 and GRB 021004 marginally violate these limits. If 
such violation is confirmed by the Swift satellite, a possible explanation is that 
(1) the compact objects that power GRB jets continue to eject an intermittent 
outflow for a very long timescale (> lday), (2) the GRB jet from the central 
engine has a temporal anisotropy with a large brightness contrast > 10 and small 
angular structure < 10~ 2 , or (3) many (> 10 3 ) regions fluctuate simultaneously 
in the emitting site. 

Subject headings: gamma rays: bursts — gamma rays: theory — relativity 



1. Introduction and summary 

The standard synchrotron shock model has been successful in explaining overall features 
of the gamma-ray burst (GRB) afterglows (e.g., Zhang & Meszaros 2004; Meszaros 2002). 
The standard afterglow model assumes a single relativistic blast wave expanding into an 
ambient medium with a spherical, smooth density distribution. The emitting surface (shock 
front) is assumed to be homogeneous and spherical. Such a model predicts smooth afterglow 
light curves. 
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However, recent dense monitoring of afterglow lightcurves indicates that at least some 
afterglows deviate from a smooth power law decay, such as in GRB 021004 (e.g., Lazzati et 
al. 2002; Kobayashi & Zhang 2003; Fox et al. 2003; Heyl & Perna 2003; Nakar, Piran, 
& Granot 2003) and GRB 030329 (e.g., Uemura et al. 2003; Lipkin et al. 2004; Torii 
et al. 2003; Sato et al. 2003; Urata et al. 2004). Also short-term variabilities (with the 
variability timescale At shorter than the observed time since the burst t) are observed in the 
afterglows of GRB 011211 (Jakobsson et al. 2004; Holland et al. 2002) and GRB 021004 
(Bersier et al. 2003; Halpern et al. 2002). These variabilities carry a wealth of information 
about the central engine and its surroundings. 

So far four major scenarios have been proposed for afterglow variabilities. (1) The 
ambient density into which a blast wave expands may have fluctuations (see § 3.1). (2) 
The emitting surface may have an intrinsic angular structure, i.e. the so called patchy shell 
model (see § 3.2). (3) The shocks may be "refreshed" by slow shells that catch up with 
the decelerated leading shell (see § 3.3). (4) The central engine may be still active at the 
observing time (see § 3.4). 

In this paper, we show that some kinds of afterglow variability are kinematically forbid- 
den under some standard assumptions. We assume that the standard model determines the 
power-law baseline of the afterglow flux F u , and derive the following limits for dips (bumps) 
that deviate below (above) the baseline with a timescale At and amplitude AF U : 

(a) No dips in afterglow light curves can have a larger amplitude than the limit given by 
equation (4) (see § 2). 

(b) Ambient density fluctuations cannot make a bump in afterglow light curves larger than 
that in equation (7) (see § 3.1). 

(c) Patchy shells cannot make a bump with a timescale shorter than the observed time 
At < t, although the rising time At rise < t is allowed (see § 3.2). 

(d) Refreshed shocks cannot make a bump with At < t/4 (see § 3.3). 

These limits are summarized in Figure 1. Notice that the limits (c) and (d) are derived from 
purely geometrical arguments, and hence, they only give constraints on At (not on AF V ). 
When many regions fluctuate simultaneously, the limits (a) and (b) are modified and are 
given by equations (Al) and (A2), respectively, as discussed in Appendix § A. 

The observed variabilities in the afterglows of GRB 011211 (Jakobsson et al. 2004; 
Holland et al. 2002) and GRB 021004 (Bersier et al. 2003; Halpern et al. 2002) may 
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actually violate all the above limits (a)-(d), although these are still within errors. We suggest 4 
that this may imply that 

(1) the central engine is still active at the observed time (> 1 day; see § 3.4), 

(2) the GRB jets have a temporal anisotropy with a small angular structure and large 
brightness contrast [equations (9) and (10) in § 4], 

(3) or many (> 10 3 ) regions fluctuate simultaneously in the emitting site (see § A). 

Therefore variabilities in afterglow light curves may provide important clues to the nature 
of the compact object that triggers the burst and its surroundings. The Swift satellite will 
significantly increase such samples, and it will allow us to further understand the nature of 
GRB engines. 



First let us consider dips in afterglow light curves at which the flux temporarily decreases 
below the expected power-law decay with a duration At smaller than the observing time t, 
i.e., At < t. We assume that dips arise from nonuniformity on the emitting surface induced 
by a certain mechanism (e.g., density fluctuations). Because of the relativistic beaming, the 
visible half-angular size of an emitting surface with a Lorentz factor 7 is about 7 -1 . Since 
the emitting surface with a radius ~ R has a curvature, two photons emitted on the line- 
of-sight and at the edge of the visible surface travel in different time to the observer. This 
angular spreading time mainly determines the observed time since the burst, i.e., t ~ R/2 / y 2 c 
(Fenimore, Madras, & Nayakshin 1996; Sari & Piran 1997). If a variable region has a half- 
angular size A9(< 7 _1 ) and a viewing angle (with respect to the center of the variable region) 
#t>(< 7 _1 ) to the observer, the angular spreading effect also puts a lower limit on the variable 
timescale 5 as At > i?A6*max[A#/2, 26 v ]/c. (This does not depend on the Lorentz factor of 
the variable region.) Then the ratio of the variable area AS* ~ ir(RA9) 2 to the visible area 
S ~ 7r(i?/7) 2 has an upper limit, 



4 Jakobsson et al. (2004) have also suggested the long-acting engine and the asymmetric density or energy 
variations as the origin of the variabilities observed in GRB 011211 through geometrical considerations. 



2. Dips in afterglow light curves 




(on — axis) 
(off — axis) 



(1) 



5 Note that At > R[(6 V + A9) 2 - (6 V - A6) 2 ]/2c = 2RA69 v /c for the off-axis case. 
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where we use the typical value 9 V ~ 7 x /2 for the off-axis case 9 V > A#. 

In addition, the thickness of the blast wave affected by a certain mechanism (e.g., density 
fluctuations) should be less than the variability timescale multiplied by the speed of light 
~ cAt. Since the overall thickness of the blast wave is about ~ -R/I67 2 ~ ct/8, the ratio of 
the variable volume AV ~ cAtAS to the visible volume V ~ ctS/ 8 is less than 

AV < r 8(At/t) 2 , (on -axis) (2) 



V ~ \ 2(At/t) 3 . (off -axis) 

To obtain the upper limits on the amplitude of a dip, we assume that the emission from the 
variable volume is suddenly shut off. Additional effects (e.g., cooling timescale) only make 
the dip less significant. Then the kinematical upper limits for the dips are 

\AF U \ f 8(At/t) 2 , (on -axis) 

~ V ~ \ 2(At/t) 3 , (off -axis) 

regardless of the cause of the variability, as long as the dips are produced by a disturbance 
on the emitting surface. 

When deriving equation (3) from equation (2), we assumed that the visible region has 
a uniform brightness. In reality a spherical afterglow appears on the sky as a ring because 
of the relativistic effect (Panaitescu & Meszaros 1998; Sari 1998; Waxman 1997). Since 
the surface brightness normalized by its mean value is about I v jl v ~ 0.1 at the center and 
about I v jl v ~ 3 at the edge in the optical band [F v oc i/ 1_ p)/ 2 ) even for a uniform surface 
(Granot & Loeb 2001), we replace equation (3) by 

1^1 < I (4/5)(At/t) 2 , (on -axis) (4) 



F v ~ \ Q(At/ty. (off -axis) 

The above limit is applicable when the emitting surface has one single variable region. 
When many regions fluctuate simultaneously, one may use equation (Al) instead, as dis- 
cussed in § A. Notice, however, that the limit in equation (4) is still useful since its violation 
implies that there are many variable regions. 



3. Bumps in afterglow light curves 

Since the surface brightness of variable regions has the lower limit (zero), we could give 
constraints on dips in a general way. However, when considering afterglow bumps, since the 
surface brightness of variable regions has no upper limit in principle, its constraint should 
depend on the specific model of bumps. In the following we separately discuss each probable 
scenario in turn. 
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3.1. Ambient density fluctuations 

Ambient density fluctuations may lead to afterglow variabilities (Wang & Loeb 2000; 
Lazzati et al. 2002; Dai & Lu 2002; Nakar, Piran, & Granot 2003). Such fluctuations 
may be due to turbulence of interstellar medium or variable winds from the progenitor star. 
Here we obtain a kinematical upper limit on the variabilities that does not depend on the 
properties of the density fluctuations. 

As in the case of dips, which we have discussed in the previous section, the ratio of the 
variable volume AV to the visible volume V satisfies equation (2). Since the kinetic energy 
in the visible volume -Ekin is almost uniformly distributed and the conversion efficiency from 
the kinetic energy E kin to the internal energy E (lab- frame) by shocks is close to unity (i.e., 
E ~ E kin ) if the shocks are relativistic, we have AE/E < AE kin /E kin ~ AV/V where 
AE(< AE k i Q ) is the internal (lab-frame) energy to produce the variability. Therefore the 
bolometric luminosity ratio of the variable part AL to the base level L is less than 

AL AE/At f 8f-\At/t), (on - axis) 

L ~ f c E/t -\2f-\At/t)\ (off -axis) [) 

where we assume L ~ e e f c E/t, which is a good approximation for the standard afterglow, 
and AL < e e AE/At, which does not depend on the precise radiative process. 6 The fraction 
of cooling energy is about f c ~ (vm/vc)^ 2 ^ 2 ~ 1/2 for the typical standard afterglow, i.e., 
the cooling frequency v c ~ 10 15 Hz, the characteristic synchrotron frequency v m ~ 10 12 Hz, 
and the electron power-law distribution index p ~ 2.2 (Sari, Piran, & Narayan 1998). 

Since the variable flux vAF v at an observed frequency v is clearly less than the bolo- 
metric variable flux AF (i.e., vAF v < AF), we can put the upper limits on bumps due to 
density fluctuations as 

< JL^L < I 8f-\F/uF u )(At/t), (on - axis) 
F v ~ vF v F ~ \ 2f- 1 (F/uF u )(At/t) 2 , (off - axis) {) 

where F is the bolometric base flux. The second inequality in equation (6) was derived by 
using AF/F < AL/L and equation (5). This is because the bolometric flux F {AF) is 
proportional to the luminosity L (AL) divided by the solid angle into which the emission 
is beamed, and the density enhancement only decelerates the emitting matter to reduce 



6 We implicitly assume that the energy fraction that goes into electrons e e is determined by the microscopic 
physics and it is constant. 
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the relativistic beaming. 7 We may estimate the factor F/i>F v in equation (6) assuming the 
standard afterglow model as FjvF v ~ (z//z/ c )^ p ~ 3 ^ 2 for v m < v < v c (the optical band at 
t ~ 1 day) and F/b>F v ~ (z//z/ c )( p ~ 2 )/ 2 for v m < v c < v (the X-ray band at £ ~ 1 day), since 
the synchrotron flux f c i<V c at the cooling frequency z/ c is about the bolometric flux F for 
p ~ 2.2 (Sari, Piran, & Narayan 1998). Since v c ~ 10 15 Hz at t ~ 1 day, we have FjvF v ~ 1 
for the optical and X-ray bands (y > 10 15 Hz). 

Finally, taking the ring-like image of the afterglow into account as in equation (4), we 
replace equation (6) by 

AF^ f (A/5)f-\F/uF u )(At/t), (on -axis) 

F„ -\Qf-\F/vF u )(At/tf. (off -axis) { ' 1 
Note that the above derivation uses only the properties of the standard afterglow. 

The above limit may be applied when the emitting site has a single variable region. 
When many regions fluctuate simultaneously, we may use equation (A2) instead, as discussed 
in § A. Note, however, that the limit in equation (7) is still useful since its violation implies 
many variable regions. 



3.2. Patchy shell model 

An intrinsic angular structure on the emitting surface (patchy shell) may also produce 
the variability of the afterglow (Meszaros, Rees, & Wijers 1998; Kumar & Piran 2000a). 
Such patchy shells may be realized in the sub-jet model (Yamazaki, Ioka, & Nakamura 2004; 
Ioka & Nakamura 2001a). Since the visible size ~ 7 _1 grows as the Lorentz factor 7 drops, 
the observed flux varies depending on the angular structure. 

The variability timescale is always At > t (Nakar & Oren 2004) if the angular fluctu- 
ation is persistent (see also § 4). The rise time of the variability At risc could be At risc /t ~ 
27 A# < 1 since it is determined by the timescale on which the angular fluctuation A9 enters 
the visible region. (The lateral expansion has to be slow for A# < 7 _1 ; see § 4.) However 
it takes ~ t for the flux to go back to its mean level since the visible region expands on the 
timescale ~ t. 



7 If we consider a void in the ambient medium, instead of the density enhancement, the matter freely 
expands in the void, so that the Lorentz factor becomes higher than that in other parts. However the 
difference within At < t is only a factor of ~ 2 and negligible for an order-of-magnitude argument. 
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3.3. Refreshed shocks 

Afterglow variabilities may also arise in the refreshed shock scenario, in which multiple 
shells are ejected instantaneously (i.e., the ejection timescale is comparable to the GRB 
duration and short compared to the observed time) but the inner shell is so slow that it 
catches up with the outer shells on a long timescale when the Lorentz factor of the outer shells 
drops slightly below that of the slow shell (Rees & Meszaros 1998; Panaitescu, Meszaros, & 
Rees 1998; Kumar & Piran 2000b; Sari & Meszaros 2000; Zhang & Meszaros 2002). Since 
inner shells only increase the afterglow energy, the observed flux does not go back to the 
original level (i.e., the extrapolation from the original power-law). 

The variability timescale should be At > t/A if the acceleration of the GRB ejecta 
is hydrodynamic. This is because, if hydrodynamically accelerated, the slow shell has an 
opening angle larger than the inverse of its Lorentz factor A6 > and a factor of ~ 2 
dispersion in the Lorentz factor of the slow shell 7 S , since the shell is hot in the acceleration 
regime and expands with a sound speed ~ c in the comoving frame. If A6 > 7" 1 , the 
variability timescale is equal to or larger than the observed time (angular spreading time) 
At > R/2cjg ~ t/A, since the Lorentz factor of the slow shell is comparable to that of the 
blast wave 7^ ~ 27 when they collide. Even when A9 < 77 1 , if the Lorentz factor of the 
slow shell 7 S varies at least by a factor of ~ 2, the slow shell spreads to have a width in the 
lab-frame A ~ R/^ll ~ ct/4, so that the variability timescale is again At > t/A. Therefore 
another acceleration mechanism is required for At < t/A. 

The slow shell may satisfy A6 < 77 1 if the outer shell has an opening angle smaller than 
77 1 . This is because a part of the slow shell outside the wake of the outer shell is decelerated 
and cut off by the ambient material. Only the shell in the wake remains cold and freely 
expands. This mechanism may explain the bumps in GRB 030329 with At ~ t j < t where 
tj is the jet break time (Granot, Nakar, & Piran 2003). However the dispersion of 7^ should 
be small for A < cAt, so that a non-hydrodynamical acceleration is still needed. 

3.4. Long-acting engine model 

A bump at an observed time t may suggest that the central engine is still active at 
that time t (Rees & Meszaros 2000; Zhang & Meszaros 2002; Dai & Lu 1998). A very long 
duration could arise if it takes a long time for the disk around a black hole to be completely 
absorbed, such as in some cases of the collapsar models (e.g., MacFadyen, Woosley, & Heger 
2001) or if the central object becomes a neutron star, such as a millisecond magnetar (e.g., 
Usov 1994). 
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Afterglow variabilities may arise when a long-term intermittent outflow from the central 
engine collides with the preceding blast wave, 8 with itself (internal shock), or with the pro- 
genitor stellar envelope. In principle, the variability timescale could be down to millisecond 
level (light crossing time of the central engine) and there is no limit on the flux variability. 
The minimum energy to produce the variability is 

AE vF v AF U At Q 

~e ~ ^"fTT^^' (8) 

where Q is the solid angle into which the variable emission is collimated and E is the afterglow 
energy in the visible region. Since the solid angle Q may be as low as ~ n% 2 , the inverse 
square of the Lorentz factor of the emitting matter, there is in principle no lower limit for Q, 
and hence for the minimum energy AE, if we consider a high Lorentz factor 7^. Considering 
other effects (such as multi-shock emission upon collision) may require a larger energy (Zhang 
& Meszaros 2002). High energy gamma-rays could be a diagnosis of the model (Ramirez- Ruiz 
2004). 



3.5. Others 

There are several other possibilities to produce variabilities in afterglow light curves: the 
gravitational microlensing (Loeb & Perna 1998; Garnavich, Loeb, & Stanek 2000; Ioka & 
Nakamura 2001b), combined reverse shock and forward shock emission (Kobayashi & Zhang 
2003), supernova bumps (e.g., Bloom et al. 1999), and dust echos (Esin & Blandford 2000; 
Meszaros& Gruzinov 2000). These are not repeated, and the variability timescales are 
usually At/t > 0.1. Since these variabilities have distinct temporal and spectral features, we 
will be able to distinguish these possibilities from the mechanisms we have discussed in this 
paper. 



4. Solutions to forbidden variabilities 

We have studied kinematical limits on the afterglow variabilities for both dips and 
bumps, which are summarized in Figure 1: (a) Dips have a smaller amplitude than that 
given by equation (4), (b) Density fluctuations cannot make a bump larger than the limit 
of equation (7), (c) Patchy shells cannot make a bump with a timescale shorter than the 



8 In the refreshed shock scenario (in the previous section), shells are ejected at almost the same time, not 
for a long time. 
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observed time At < t, (d) Refreshed shocks cannot make a bump with At < t/4. If many- 
regions fluctuate simultaneously, the limits (a) and (b) are replaced by equations (Al) and 
(A2), respectively, as discussed in § A. The variabilities in GRB 011211 (At/t ~ 0.1 and 
\AF V \/F V ~ 0.1; Jakobsson et al. 2004; Holland et al. 2002) may violate the limits (a), (c) 
and (d) and marginally violate (b), while the variabilities in GRB 021004 (At/t ~ 0.01 and 
\AF U \/F U ~ 0.05; Bersier et al. 2003; Halpern et al. 2002) may violate all these limits, if a 
single region fluctuates, although still within errors. 

(1) One possible explanation to the forbidden variabilities is the day-long central engine 
model discussed in § 3.4. Interestingly, this scenario may easily produce metal features (Rees 
& Meszaros 2000) and the metal emission lines are indeed observed in the X-ray afterglow 
of GRB 011211 (Reeves et al. 2002). A simple form of this model may explain forbidden 
bumps, but not dips. Nevertheless, as noted by Rees & Meszaros (2000), we cannot rule out 
the possibility that all afterglow is due to the central engine itself. This extreme version of 
this model may resolve the forbidden bumps as well as dips in the afterglow. 

(2) Another solution could be provided by non-standard assumptions that the emitting 
surface is anisotropic (during the period when light curves smoothly decay) and that the 
anisotropy is temporal. 9 Since equation (1) is derived from geometrical arguments and 
applicable to both dips and bumps, the angular size of the anisotropy needs to be 

Afl< J T^t/ty/*- 10- 2 , (on -axis) 

- \ (l/2) 7 - 1 (At/t) ~ 10- 3 , (off - axis) 1 ' 

and with the temporal timescale ~ At the surface brightness in this region should be en- 
hanced (bumps) or reduced (dips) by 

\AI V \ _ \AF V \ S J (|AF„|/F„)(t/At) ~ 10, (on - axis) 

— 1 /I/I A 77 1 I I \/-#./A-#.\2 1A3 /„« -„\ V 1U / 



F v AS ~ 1 A{\AF V \/F v ){t/Atf ~ 10 3 , (off - axis) 

where numerical values are for GRB 021004 (At/t ~ 0.01 and \AF U \/F U ~ 0.05) and we have 
assumed 7 ~ 10 at t ~ 1 day. For the off-axis case, the brightness contrast may be too large 
\AI V \/I U ~ 10 3 to be reconciled with the observed narrow distribution of the geometrically 
corrected gamma-ray energy (Frail et al. 2001). 

Let us examine each violation one by one in this temporally anisotropic model. To 
explain narrow dips violating (a), the variable region should be initially brighter than the 
limit given by equation (10) and then it should be darkened, such as by the density fluctu- 
ations. However the angular size should also satisfy equation (9). This size is contrary to 



9 In the patchy shell model, the anisotropy is not temporal but persistent (see § 3.2). 
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our common belief that the initial fluctuations with A9 < 7"" 1 are erased during the fireball 
evolution because the visible region ~ 7 _1 is causally connected. The lateral expansion has 
to be slower than the usual assumption (i.e., the sound speed in the local frame) for the 
initial fluctuations to prevail. [Notice that several numerical simulations imply a slow lateral 
expansion (Granot et al. 2001; Kumar & Granot 2003; Cannizzo, Gehrels, & Vishniac 
2004).] We also expect that other variabilities due to the patchy shell effect should appear 
in the light curves when the angular fluctuations enter the visible region. 

Next let us consider bumps violating the limits (b)-(d) in the temporally anisotropic 
model. To explain variabilities violating the limit (b) in the density fluctuation scenario, the 
angular energy distribution has to be initially anisotropic and the energetic spot needs to be 
brightened by the density fluctuation, as in the dip case (a). Again the lateral expansion has 
to be slow, and we also expect other variabilities due to the patchy shell effect. In the patchy 
shell model violating (c), the anisotropy should be temporal ~ At, so that an external factor 
such as density fluctuations may be necessary. Also in this case the lateral expansion has 
to be slow. In order for the bumps due to refreshed shocks to violate (d), the acceleration 
mechanism should be non-hydrodynamical (see § 3.3). If so, the rise time At r ise may be as 
short as At rise /t ~ (7A#) 2 < 1 (on-axis), but the flux does not return to the original level. 
The step-wise light curve is a signature of the refreshed shock. 

(3) The other solution could be that many regions fluctuate simultaneously. As discussed 
in § A, the limits (a) and (b) are modified in this case and are given by equations (Al) and 
(A2), respectively. The limit (a) for dips in equation (Al) may be still violated by GRB 
011211 and GRB 021004, while the limit (b) for bumps due to density fluctuations in equation 
(A2) may not be violated. Even so the number of variable regions has to be larger than ~ 10 3 
to reconcile the limit (b) with GRB 021004, and this suggests that the mean separation of 
the density clumps (with a radius ~ 10 14 cm) is about ~ 10 15 cm (see § A). Therefore these 
limits provide interesting constraints on density fluctuations. In this model the anisotropy of 
the emitting surface is also strong and temporal as in the previous model, and the observed 
flux almost always differs from the base level. 

(4) The violations could be attributed to uncertainties of the observations and data 
analyses. The violation of equation (4), which gives constraints on dips, may be due to the 
fitting scheme, because spurious dips might be produced by inappropriate power-law fitting 
to a light curve containing bumps. Intensive afterglow monitorings such as those by the Swift 
satellite and its ground-based follow-up observations, when combined with appropriate data 
fitting methods, might be able to lead to a better determination of the base level of the 
decay, which could be used to verify or refute the presence of the forbidden variabilities. 

Some afterglows show little variabilities despite dense sampling (e.g., Laursen & Stanek 
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2003; Gorosabel et al. 2004; Stanek et al. 1999). The variety could arise from the viewing 
angle (e.g., the anisotropy depends on the viewing angle) or the intrinsic diversity (e.g., each 
burst has a different anisotropy), but future observations are needed to fix its origin. 

We have considered a fixed lab-frame time to relate the emitting area and volume with 
the observed time in equations (1) and (2). If we take the equal arrival time surface into 
account (Panaitescu & Meszaros 1998; Sari 1998; Waxman 1997), these relations would be 
modified and make the limits more stringent especially near the limb of the afterglow image. 
This is also an interesting future problem. 

We thank P. Meszaros, D. Lazzati and the referee for useful comments. This work 
was supported in part by the Eberly Research Funds of Penn State and by the Center for 
Gravitational Wave Physics under grants PHY-01-14375 (for KI, SK), NASA Swift Cycle 1 
GI Program (for SK, BZ) and NASA NNG04GD51G (for BZ). 

A. The case of many variable regions 

So far we have considered only a single variable region to make the discussions simple. In 
reality many regions could fluctuate simultaneously (e.g., Gruzinov & Waxman 1999). This 
leads to weaker limits for dips (a) and bumps due to density fluctuations (b) in equations 
(4) and (7), respectively. In any case, we can still obtain meaningful limits by extending the 
previous arguments, as shown below. 

Let N v be the mean number of variable regions contributing in the observed time interval 
of ~ At. Then the flux deviates from the baseline ~ N v times larger than the single variable 
region case. However the flux variability is determined by the differences between time 
intervals. This is determined by the variance of the number of variable regions from one 
time interval ~ At to the next, and it is ~ N l J 2 if we assume the Poisson statistics. The 
Poisson statistics may be applied since different variable regions fluctuate at different lab- 
frame times. Therefore the limits for a single variable region in equations (4) and (7) have 

1/2 

to be multiplied by ~ N v ' when many regions fluctuate. 

A.l. Limits (a): dips 

Let us consider the maximum of the number of variable regions N v in the dip case (a). 
In the on-axis case, the maximum number is about the ratio of the overall thickness of the 
blast wave ~ -R/I67 2 ~ ct/S to the fluctuating thickness ~ cAt, i.e., N v < (t/At)/S. In 



the off-axis case, it is about the ratio of the visible volume V to the variable volume Ay in 
equation (2), i.e., N v < V/AV ~ (t/Atf/2. 

Then, multiplying the limits for a single variable region in equation (4) by ~ N^ 2 , we 
obtain the limits for dips due to many variable regions as 

|AF„| f ( v / 2/5)(At/t) 3 / 2 , (on -axis) 
F v ~ \ (6/ v / 2)(At/t) 3 / 2 . (off -axis) 

The off-axis limit is shown by the dashed line in Figure 1. The variabilities in GRB 011211 
(At/t ~ 0.1 and \AF U \/F U ~ 0.1) and GRB 021004 (At ft ~ 0.01 and \AF U \/F U ~ 0.05) may 
still violate these limits. 

The above limits are quite robust, but may be too strict because the interior of the blast 
wave also needs to fluctuate in the same pace as the front and back of the blast wave when 
the number of variable regions N v is nearly maximum. A reasonable fluctuating site may 
be the front and back of the blast wave that are affected by the density clumps or the inner 
shells. If this is the case, the maximum of N v is ~ 1 in the on-axis case, and is about the ratio 
of the visible area S to the variable area AS in equation (1), i.e., N v < S/AS ~ 4(t/At) 2 , 
in the off-axis case. 



(Al) 



A. 2. Limits (b): bumps due to ambient density fluctuations 

When the density fluctuations make bumps in the afterglow, the bumps are mainly 
produced at the shock front (i.e., not such as in the interior of the blast wave). Then the 
maximum of the number of variable regions N v is ~ 1 in the on-axis case. In the off-axis 
case, it is about the ratio of the visible area S to the variable area AS in equation (1), i.e., 
N v < S/AS ~ A(t/At) 2 . 

Multiplying the limits for a single variable region in equation (7) by ~ N^ 2 , we obtain 
the limits for bumps due to many density fluctuations as 

f (A/5)f-\F/uF v )(At/t), (on-axis) 
F v ~ \ 12 f-\F/vF v ){At/t). (off -axis) 1 ' 

The off-axis limit is shown by the dashed line in Figure 1. This maximum off-axis limit is 
not violated by the variabilities in GRB 011211 (At/t ~ 0.1 and \AF U \/F„ ~ 0.1) and GRB 
021004 (At/t ~ 0.01 and \AF U \/F U ~ 0.05). Still we need N v > 1600 variable regions to 
reconcile the off-axis limit with GRB 021004 because the variabilities in GRB 021004 are 
~ 40 ~ V4600 times larger than the off- axis limit for a single variable region in equation (7). 
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This corresponds to the mean separation of the density clumps ~ [7r(ii ) /7) 2 cAt7 2 /iVj 1 / 3 ~ 
10 15 cm since the shock front with a radius R/^y ~ ct^y ~ 10 16 cm sweeps a distance cAt'j 2 ~ 
10 15 cm for t ~ 1 day, At ~ 10 3 sec and 7 ~ 10. Note that we can also estimate the clump 
radius as cAt7 ~ 10 14 cm. In this way, the limits in equations (7) and (A2) can impose 
interesting constraints on the density fluctuations. 
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Fig. 1. — Kinematically allowed regions for afterglow variabilities are shown in the plane of 
the relative variability timescale At/t and relative variability flux AF U /F U . We have four 
limits: (a) equation (4) for dips, (b) equation (7) for bumps due to density fluctuations, (c) 
At >t for bumps due to patchy shells and (d) At > t/4 for bumps due to refreshed shocks. 
For (a) and (b), the on-axis cases are shown. When many regions fluctuate simultaneously, 
the limits (a) and (b) are replaced by equations (Al) and (A2), respectively, and the off-axis 
cases are shown by dashed lines. We assume FjvF v ~ 1 and f c ~ 1/2. The variabilities 
in the afterglows of GRB 011211 (Jakobsson et al. 2004; Holland et al. 2002) and GRB 
021004 (Bersier et al. 2003; Halpern et al. 2002) are also shown by dotted circles, and may 
violate some of these limits. 



